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RESEARCHES IN THE LUNAR THEORY. 
By Gr. W. Hill, Nyack Turnpike, N. Y. 

(Continued from p. 26.) 

Chapter II. 

Determination of the inequalities which, depend only on the ratio of the mean motions 

of the sun and moon. 



If the path of a body, whose motion satisfies the differential equations 
d?x 
di 2 



— 2n' \ + [t — 3n' 2 ~\ x = 0, 
dt i-r J 

§+2»'~+ y =0, 

dt' dt r 

intersect the axis of x at right angles, the circumstances of motion, before 
and after the intersection, are identical, but in reverse order with respect to 
the time. That is, if t be counted from the epoch when the body is on the 
axis of x, we shall have 

x — function (t 2 ) , y — t . function (t 2 ) . 

For if, in the differential equations, the signs of y and t are reversed, but that 
of x left unchanged, the equations are the same as at first. 

A similar thing is true if the path intersect the axis of y at right angles ; 
for if the signs of x and t are reversed, while that of y is not altered, the 
equations undergo no change. 

Now it is evident that the body may start from a given point on, and at 
right angles to, the axis of x, with different velocities ; and that, within 
certain limits, it may reach the axis of y, and cross the same at correspond- 
ingly different angles. If the right angle lie between some of these, we 
judge, from the principle of continuity, that there is some intermediate 
velocity with which the body would arrive at and cross the axis of y at right 



angles. 
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The difficulty of this question does not permit its being treated by a 
literal analysis ; but the tracing of the path of the body, in numerous 
special cases, by the application of mechanical quadratures to the differential 
equations, enables us to state the following circumstances : — 

If the body be projected at right angles to, and from a point on, the 

the axis of x; whose distance from the origin is less than 0.33 . . . \-^ , there 

is at least one (near the limit there are two) value of the initial velocity, with 
which the body, in arriving at the axis of y, will cross it at right angles. 
Beyond this limit it appears no initial velocity will serve to make the body 
reach the axis of y under the stated condition. 

If the body move from one axis to the other and cross both of them per- 
pendicularly, it is plain, from the preceding developments, that its orbit will 
be a closed curve symmetrical with respect to both axes. Thus is obtained a 
particular solution of the differential equations. While the general integrals 
involve four arbitrary constants, this solution, it is plain, has but two, which 
may be taken to be the distance from the origin at which the body crosses the 
axis of x and the time of crossing. 

Certain considerations, connected with the employment of Fourier's 
Theorem and the po?sibility of developing functions in infinite series of 
periodic terms, show that, in this solution, the co-ordinates of the body can 
be represented, in a convergent manner, by series of the following form, 
x = A cos [y (t — £„)] + A t cos B[y(t — £„)] + A 2 cos 5 \y (t — $>)] + . • • , 
y = B sin [y (t — £„)] + B x sin 3 [y (t — $>)] + B 2 sin 5 j> (t — $,)] + . . . , 

where t denotes the time the body crosses the axis of x, and — is the time 

V 

of a complete revolution of the body about the origin. We may regard v 
and to as the arbitrary constants introduced by integration ; the coefficients 
A , Ai . . . Bi, B . . . are functions of ^, n' and v. 
For convenience sake we may put 

Ai = a, : -\- a_^_! , Bi — a, t — a_i_i . 

Then, t being put for v (t — t ), the series, given above, may be written 

x = 2*. K cos (2i + 1) t, 
y = Xi ■ &i sin (2i + 1) t , 
the summation being extended to all integral values positive and negative, 
zero included, for i. By adopting polar co-ordinates such that 

x = r cos $ , y = r sin 4> , 
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and writing v for <p — t, that is, for the excess of the true over the mean 
longitude of the moon, the last equations are equivalent to 

r cos v = X t . &i cos 2vv, 
r sin v = 2«. a^ sin 2it. 

In order to avoid the multiplication of series of sines and cosines, and 
reduce everything to an algebraic form, for x and y, we substitute the imagi- 
nary variables u and s, and put £ = e 1 v ~ 1 . We have then 

£ will always be employed as the independent variable in place of t or r. 

Denoting the operation £ — = — V — 1— by the symbol D, so that, in 

general, 

D «) = i<> 
and taking the liberty of separating this symbol as if it were a multiplier, 

and moreover putting 

n' n' u 

m = - = ■ , x = £ , 

the differential equations, determining u and s, given in the preceding chap- 
ter, may be written 

|> + 2mD + §m 2 - * 1 u + §m 2 * = 0, 

L 2 (ms)*J 2 

\l) 2 — 2mD + f m 2 — ---1 s + f m 2 « = . 
L 2 (us)H 2 

It will be noticed that either of these equations can be derived from the 
other by interchanging u and s and reversing the sign of m or D. We may 
also remind the reader that they determine rigorously all the parts of the 
lunar coordinates which depend only on the ratio of the mean motions of the 
sun and moon and on the lunar eccentricity. The Jacobian integral, in the 
present notation, is 

Du.Vs+ f%- + 4- m 2 (u + s) 2 = C. 

(us) i 4 

The most ready method of getting the values of the coefficients a {} is that 
of undetermined coefficients; the values of u and s, expressed by the pre- 
ceding summations with reference to i, being substituted in the differential 
equations, the resulting coefficient of each power of £, in the left members, is 
equated to zero, which furnishes a series of equations of condition sufficient 
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to determine all the quantities a 4 . For this purpose we may evidently employ 
any two independent combinations of the three equations last written, and 
it will be advisable to form these combinations in such a manner that the 
process of deriving the equations of condition may be facilitated in the 
largest degree. Now it will be recognized that the presence of the term 

•7-^— , in one of the factors of the differential equations, is a hindrance to 
(us)§ 

their ready integration, being the single thing which prevents them from 

being linear with constant coefficients. Hence we avail ourselves of the 

possibility of eliminating it. Multiplying the first differential equation by s, 

and the second by u, and taking, in succession, the sum and difference, 

uB 2 s + sD 2 u — 2m (uDs — sl)u) — ~ + -| m 2 {u + s) 2 = 0, 

(us)i 2 

uD 2 s — sl) 2 u — 2m (uDs + sDu) + % m 2 (u 2 — s 2 ) = 0, 

then, adding to the first of these the integral equation, and retaining the 
second as it is, we have, as the final differential equations to be employed, 

D 2 (us) —Du.Ds —2m {uDs — sDu) + -| m 2 (u + s) 2 = C, 

I) (uDs — sDu — 2mus) + ^-m 2 (u 2 — s 2 ) =0. 

2 

It must be pointed out, however, that these equations are not, in all 
respects, a complete substitute for the original equations. It will be seen that 
ft or x, an essential element in the problem, has disappeared from them, and 
that, in integration, an arbitrary constant, in excess of those admissible, will 
present itself. This will be eliminated by substituting the integrals found in 
one of the original differential equations, in which n or % is present ; the 
result being an equation of condition by which the superfluous constant can 
be expressed in terms of ft and the remaining constants. 

We remark that the left members of our differential equations are 
homogeneous and of two dimensions with respect to u and s. If the first 
were differentiated, the constant G would disappear, and both equations would 
be homogeneous in all their terms. This property renders them exceedingly 
useful when equations of condition are to be obtained between the coefficients 
of the different periodic terms of the lunar co-ordinates, and it is for this 
purpose that we have given them their present form. 
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From the signification of the symbol D, 

Du = X i . {2% + 1) a,^ + 1 , Ds - 2, . (2i + 1) a_^ 2i+1 , 

D 2 W = 2* . (2i + l) 2 a^ 21 + 1 , V 2 s = 2 i . (2i + l) 2 a_ 4 _^ 2i + * ; 

also 

us = 2,- . [2* . a^a-i.y] £ 2j ' , 
# 2 = 2, • [2,: . a i a_ i +j _ 1 ~] % y , 
s —Zj. [2* . a i a_ i _y_ 1 J c, 3 , 
Du.Ds = — 2, . [2, , (2» + 1) (2» - 2? + 1) a,a,_ J £* , 
w 2) s _ S X> W - _ 22,- . [2, . (2i — j + 1) a,a,_J ^ , 

where the summations with reference to j have the same extension as those 
with reference to i. On substituting these expressions in the differential 
equations, and equating the general coefficients of Z^ to zero, we get 

2,. [(2* + 1) (2i - 2j + 1) + 4f + 4 (2» — j + 1) m + -| m 2 ] a,a w 
+ -- m 2 2*. [a i a_ i+ ^_ 1 + a i a_ j _^_ 1 J = 0, 

4j X t . \2i — j + 1 + mj a«a w — -=- m 2 2«. [a i a_ i+ ^_ 1 — a i a_ i _ / _ 1 J = 0, 

which hold for all integral values of j both positive and negative except that, 
when j — 0, the right member of the first equation is C instead of ; but as 
the second equation is an identity for j ■= 0, for the present this value of j 
will be excluded from consideration. 

By multiplying the first equation by 2, and the second by 3, and taking 
in succession the difference and sum, the simpler forms are obtained, 

2,. [8i 2 — 8 (4j — 1) i + 20j 2 — 16? + 2 + 4 (4» — 5j + 2) m + 9m 2 ] a A _ y 

+ 9m 2 2*.a i a_ i+i _ 1 = 0, 
2,. [8i 2 + 8 (2j + l)i - 4f + 8j + 2 + 4(4i + j + 2) m + 9m 2 ] a^--,- 

+ 9m 2 2«.a i a_ i _^_ 1 = 0. 

These two equations are not distinct from each other, when negative, as well 
as positive values, are attributed to j. For if, in the expression under the first 
sign of summation in the first equation, we substitute, which is allowable, for 
j ( i — j t and — j forj throughout the equation, the result is identical with the 
second equation. This is explained by the fact that we get all the inde- 
pendent equations of condition, these equations are capable of furnishing, by 
attributing only positive values to j. Hence, allowing j to receive positive 
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and negative values, all the equations of condition can be represented by a 
unique formula. 

Although the number of these equations is infinite, and also that of the 
coefficients a, { , it is not difficult to see that the first ought to be regarded as 
one less than the second; and that, in consequence of the bi-dimensional 
character of the equations, they suffice to determine the ratio of any two of 
the quantities ^ in terms of m. It will be seen, from developments to be 
given shortly, that if m is regarded as a small quantity of the first order, a 4 is 
of the ± 2i ih order. It will be advisable then to select a as the coefficient to 
which to refer all the rest ; and we shall have, in general, 

a; — a F(m). 
The equations of condition, as written above, determine the a, in pairs ; that 
is, if we put j = 1, we have the equations suitable for determining a x and a_ x , 
and, in general, the equations, as written, determine a,- and a_,- . And, as 
they involve both these quantities, it will be advantageous to eliminate 
approximately each in succession, as far as that can be done without depriv- 
ing the equations of their bi-dimensional character. 

By putting, in succession, in the terms under the first sign of summation, 
i = and i = j, it will be found that these equations contain, severally, the 
terms 
[20/ — 16j + 2 — 4 (5j — 2) m + 9m 2 ] a a_,. 

+ [- 4f- 8j + 2 - 4 (j-2) m + 9m 2 ] a,a, , 
[— 4f + Sj + 2 + 4 (j + 2) m + 9m 2 ] a a_, 

+ [20f + 16/ + 2 + 4(5; + 2)m + 9m 2 ] a a, , 
which are the terms of principal moment in determining a_,- and a,- . Let us 
then multiply the first equation by 

— 4/ + 8J + 2 + 4(j + 2) m+ 9m 2 , 
and the second by 

— 20f + 16j — 2 + 4 (5j — 2) m — 9m 2 , 

and, adding the products, divide the whole by 

48f [2 (4f — 1) — 4m + m 2 ] . 
Then, adopting the notation 

i 4 {j — 1) i + 4f + 4j — 2 — 4 (i — j + 1) m + m 2 



1***1 



j 2 (4f — 1) — 4m + ni 



2 



f 1 3m^ 4f — 8j — 2 — 4 ( j + 2) m — 9m 2 

L,?J ~~ 16f 2 (4f — 1) — 4m + m 2 
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n - — 20/ — 16j + 2 — 4 (5j — 2) m + 9m 2 
Uj " 16j 2 2 (4f — 1) — 4m + m 2 ~ ' ' 

the system of equations, which determines the coefficients a ; , is represented 
by the unique formula 

2*- [D"> *l Wi-i + U] ^a-i+j-i + (j) a 1 a^ < _ y _ 1 ] = 0, 
where j must receive negative as well as positive values. It will be perceived 
that 

[j,o] = o, Ij.il =-i; 

hence the last equation is in a form suitable for determining the value of a,-. 
The quantities \J,t], [j] and (j) admit of being expressed in a simpler 
manner ; thus 

[j {i^_i + 4i U — i) j—l — m 

U ' J j j ~ 2 (4f — 1) — 4m + m 2 ' 

whence 

LMJ + L ^ »J- ^ +2(^-1) -4m + m 2 ' 

U»*J L ^ *J- . 2 (4f — 1) — 4m + m 2 ' 

in addition 

[j] = ^. m >_3 lQf— 2j— 5 — (i+H)m ml 



16f 4j 2 2 (4/ — I) — 4m + m 2 

(j) - ... 27 m , , 3 13f+4j-5 + (5j-ll)m 
KJ} 16/ ^4f 2(4f — 1) — 4m + m 2 



m 2 , 



rfi + (_,-) - _i 3i + i+ 2m m * 

L/J tl ii- 2 j 2 ( ^ 2 _ 1) _ 4m + m2 I 

rfl _ f_fl _ 27 m2 _ £ 16/-3i-5-(3j+ll)m m , 
m v J) 8f 2f 2 (4f — 1) — 4m + m 2 

In making a first approximation to the values of the coefficients, one of 
the terms of the equation may be omitted ; for, when j is positive, the term 
2j. (j) a i a_ i _ J _ 1 is a quantity four orders higher than that of the terms of the 
lowest order contained in the equation ; and, when j is negative, the same 
thing is true of £ 4 . [j] a i a_ i+y _ 1 . Hence, with this limitation, the equation 
may be written in the two forms 

2i- [D'> CI a A-,- + Tj] a i a_ i+ ,-_ 1 ] = , 
2i- [[— j, *] a f a i+ i + (— i) aia.i+j.t] = . 
where j takes only positive values. 
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From these two equations, by omitting all terms but those of the lowest 
order, we derive the following series of equations, determining the coefficients 
to the first degree of approximation, 
a a t = [1] a a , 
a a_ 1 = ( — 1) a a , 

a a 2 = [2] [a a x + a^] + [2, 1] a x a_! , 
a a_ 2 = (— 2) [a a t + a^o] + [— 2,— 1] a^j , 
a a 3 = [3] [a a 2 + a^ + a 2 a ] + [3, 1] a x a_ 2 + [3, 2] a 2 a_j , 
a a_ 3 = ( — 3)[a a 2 + a^ + a 2 a ] + [ — 3, — 1] a_ x a 2 + [ — 3, — 2] a_ 2 a! , 
a a 4 = [4] [a a 8 + a x a 2 + a^ +a 3 a ] + [4, 1] a^g + [4, 2] a 2 a_ 2 + [4, 3] a 3 a_! , 
a a_ 4 = (—4) [a a 3 + a t a 2 + a 2 a x + a 3 a ] -f [ — 4, — 1] a^a,, + [— 4, — 2] a_ 2 a 2 
+ [— 4, — 3]a_ 3 a 1 , 

The law of these equations is quite apparent, and they can easily be extended 
as far as desired. The first two give the values of a x and a_ lt the following 
two the .values of a 2 and a_ 2 by means of the values of a x and a_ x already 
obtained, and so on, every two equations of the series giving the values of 
two coefficients by means of the values of all those which precede in the order 
of enumeration. A glance at the composition of these equations must con- 
vince us that all attempts to write explicitly, even this approximate value 
of n-i, would be unsuccessful on account of the excessive multiplicity of the 
terms. However, they may be regarded, in some sense, "as giving the law of 
this approximate solution, since they exhibit clearly the mode in which each 
coefficient depends on all those which precede it. As to the degree of approxi- 
mation afforded by these equations, when the values are expanded in series 
of ascending powers of m, the first four terms are obtained correctly in the 
case of each coefficient. Thus a x and a_ x are affected with errors of the 6th 
order, a 2 and a_ 2 with errors of the 8th order, a 3 and a_ 3 with errors of the 
10th order, and so on. 

The values of these quantities can be determined either in the literal 
form, where the parameter m is left indeterminate, as has been done by Plana 
and Delaunay, or as numbers, which mode has been followed by all the 
earlier lunar-theorists and Hansen. In the latter case, one will begin by 
computing the numerical values of the quantities [j, i], [j] and (j), corre- 
sponding to the assumed value of m. for all necessary values of the integers 
i and j. 
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The great advantage of our equations consists in this, that we are able 
to extend the approximation as far as we wish, simply by writing explicitly 
the terms, our symbols giving the law of the coefficients. How rapid is the 
approximation in the terms of these equations will be apparent, when we say, 
that, after a certain number of terms are written, in order to carry this four 
orders higher, it is necessary to add to each of them only four new terms ; 
and, thereafter, every addition of four terms enables us to carry the approxi- 
mation four orders farther. 

The process which may be followed to obtain the values of the a 4 with 
any desired degree of accuracy, is this: — the first approximate values will be 
got from the preceding group of equations until the a^ become of orders 
intended to be neglected ; then one will recommence at the beginning, using 
the equations each augmented by the terms necessary to carry the approxi- 
mation four orders higher; substituting in the new terms the values obtained 
from the first approximation, and, in the old, ascertaining what changes are 
produced by employing the more exact values instead of the first approxi- 
mations. A second return to the beginning of the work will, in like manner, 
push the degree of exactitude four orders higher. In this way any required 
degree of approximation may be attained. 

Whatever advantage the present process may have over those previously 
employed is plainly due to the use of the indeterminate integers i and j, 
which, although much used in the planetary theories, no one seems to have 
thought of introducing into the lunar theory. This enables us to perform a 
large mass of operations once for all. 

For the purpose of making evident the preceding assertions, and because 
we shall have occasion to use them, we write below the equations determining 
the coefficients a { correct to quantities of the 13th order inclusive. 

a ai = [1] [a? + 2a l _ 1 a 1 + 2a_ 2 a 2 ] + (1) [ai x + 2a a_ 2 + 2a!a_ 3 ] 

+ [1, — 2] a_ 2 a_ 3 + [1, — 1] a_ia_ 2 + [1, 2] a^ + [1, 3] a 3 a 2 , 

a a_ 1 = [— 1] [alt + 2a a_ 2 + 2a!a_ 8 ] + (— 1) [ag + 2a_ x a! + 2a_ 2 a 2 ] 

+ [— 1,— 3]a_ 3 a_ 2 +[— 1,— 2]a_ 2 a_!+ [—1,1] aA-f [—1, 2] a 2 a 3 , 
a a 2 = [2] [2a a x + 2a_ x a 2 + 2a_ 2 a 3 ] + (2) [2a_ x a_ 2 + 2a a_ 3 + 2a 1 a,_ i ~] 

+ [2,— 2]a_ 2 a_ 4 + [2— l]a_ ia _ 3 + [2,1] a^-f [2,3]a 8 a 1 + [2,4] a 4 a 2 , 
a a_ 2 = [ — 2] [2a_ja_ 2 + 2a a_ 3 + 2a x a_ 4 ] + ( — 2) [2a aj + 2a_ x a 2 + 2a_ 2 a 3 ] 

+ [—2,-4] a_4a_ 2 + [— 2— 3]a_ 3 a_ 1 + [— 2,— 1] a_ iai + [—2, 1] a x a 3 

+ [— 2, 2] a 2 a 4 , 
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a a 3 = [3] [a 2 + 2a a 2 -f 2a_ 1 a 3 ] + (3) [al 2 + 2a_!a_ 3 + 2a a_ 4 ] 

+ [ 3 > — !] a_ja_ 4 + [3, 1] a x a_ 2 + [3, 2] a 2 a_j + [3, 4] a^ , 

a a_ 3 = L— 3] [a!_ 2 + 2a_ 1 a_ 3 + 2a a_ 4 ] + (— 3) [a 2 + 2a a 2 + 2a_ 1 a 3 ] 

+ [ — 3,-4] a_ 4 a_!+ [ — 3, — 2] a_ 2 a x + [—3, — 1] a_ ] a 2 + [ — 3, 1] a t a 4 , 

a a 4 = [4] [2aja 2 + 2a a 3 + 2a_ x a 4 ] + (4) [2a__ 2 a_ 3 + 2a_ia_ 4 + 2a a_ 6 ] 

+ [4,-1] a_ 1 a_ 5 + [4, 1] a x a_ 3 + [4, 2] a 2 a_ 2 + [4,3] a 3 a_ x + [4,5] a^, 

a a_ 4 = [ — 4] [2a_ 2 a_ 8 + 2a_ t a_ 4 + 2a a_ 5 ] + ( — 4) [2aja 2 + 2a a 3 + 2a_ t a 4 ] 

+ [—4,-5] a_5a_H- [—4,-3] a_ 8 a x + [—4,-2] a_ 2 a 2 + [—4,-1] a^a, 

+ [—4, 1] a,a, , 

a a 5 = [5] [al + 2a.!a 3 + 2a, a 4 ] 

+ [5, 1] aja_ 4 + [5, 2] a 2 a_ 3 + [5, 3] a 3 a_ 2 + [5, 4] a 4 a_j , 

a a_ 6 = ( — 5) [al + 2a!a 3 + 2a a 4 ] 

+ [—5,-4] a_ 4 a t + [—5,-3] a_ga 2 + [—5,-2] a_ 2 a 3 + [—5,-1] a_ t a 4 , 

a a 6 = [6] [2a 2 a 3 + 2m* + 2a a 5 ] 

+ [6, 1] *&_ s + [6,2] a 2 a_ 4 + [6, 3] a 3 a_ 3 + [6, 4] a 4 a_ 2 + [6, 5] a 5 a_! , 

a a 6 — (— 6) [2a 2 a 3 + 2a,a 4 + 2a a 5 ] 

+ [_6,_5] a_ e aH- [—6,-4] a_ 4 a 2 + [—6,-3] a_ 3 a 3 + [—6,-2] a_ 2 a 4 
_|_[_e,_l]a_ 1 a 6 . 

In the first approximation 

ai = [1] a , 
a_!= (— l)a , 

a 2 =[l][2(2) + [2,l](-l)]a , 

a_ 2 = [1] [2 (-2) + [-2,-1] (-1)] a , 

or, explicitly in terms of m, 



_ 3 6 + 12m + 9m _, 



2 

nra , 



1 - 16 6 — 4m + m 2 

3 38 + 28m + 9m 2 2 

a_ t — — —^ — -~ -. ; 5— m a , 

16 6 — 4m -f- m 2 

and, after some reductions, 

27 2 + 4m + 3m 2 f^-^m-J-Qm 2 op 29 — 35m 1 4 

a2 ^^[6 -4m4-m 2 ][30-4m+m 2 ] L 238+40m+9m ~ 32 6^4lnl^J mao ' 



-27 2 + 4m + 3m 2 r 2 8-7m + 24 7 ~ m 1, 

- 2 ~64 [6 — 4m + m 2 ] [30 — 4m + m 2 ] L ^ 6— 4m + m 2 J 
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It is evident that, however far the approximation may be carried, the 
only quantities, involved as divisors in the values of the a*, are the trinomials, 
whose general expression is 

2 (4f — 1) — 4m + m 2 , 
or, particularizing, the series of divisors is 

6 — 4m + m 2 , 
30 — 4m + m 2 , 
70 — 4m + m 2 , 



It will be remarked that they differ only in their first terms, which are inde- 
pendent of m. Hence any expression, involving several divisors, can always 
be separated into several parts, each involving only one divisor, without any 
actual division by a trinomial in m. For instance, 

1 _. 1^ 1 JL_ 1 

[6 _ 4m + m 2 ] [30 — 4m + m 2 ] ~ 24 6 — 4m + m 2 24 30 — 4m + m 2 ' 

1 _ 1 1 1 _JL 

[6 _ 4m + m 2 ] 2 [30 — 4m + m 2 ] ~ 24 [6 — 4m + m 2 ] 2 24 2 6 — 4m + m 2 

H - I 1 

^ 24 2 30 — 4m + m 2 ' 

Moreover when, after this transformation, any numerator contains more or 
other powers of m than two consecutive powers, it is clear it may be reduced 
so as to contain only these by eliminating the higher powers through sub- 
tracting certain multiples of the divisor which appears in the denominator, 
or, in other words, the fraction may be treated as if it were improper. 
From this we gather that the value of a t can be expressed thus 
a< _ M , My M 2 M s 

d 

■A fll h ±1* I ^ + 

1 on ,1 — l ™2 • TOA A 1 212 I rOA A^ I 2T3 I 



6 — 


-4m -(- m 2 


30- 


-4m + m 2 
Pi 



[6- 


-4m + 


m 2 ] 2 


[30- 


-4m + 
P. 


m 2 ] 2 



[6- 


- 4m + m 2 ] 3 


[30- 


- 4m + m 2 ] 3 
P* 



+ 70— 4m + m 2 + [70 — 4m + m 2 ] 2 + [70— 4m + m 2 ] 3 + ' " * 

where M , M t . . . JVj, N 2 . . P 1} P 2 . . . are entire functions of m each of the 
form lm k + Bm k + 1 . 

The advantage of this method of treatment consists in that nothing, which is 
given by the successive approximations, would be lost, as must be the case 
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when the values are expanded in series of ascending powers of m. The 
preceding expressions, when put into this form, become 
a 1 + a_ 1 _ 2 + m 2 

— O — III , 

a 6 — 4m + m 2 

ai — a,, _ o T 9 4 — 7m ~\ 

a -°Lq 6 — 4m + m 2 J ' 

a 2 + a_ 2 _ 3 [~243 323 + 109m 9Q 23 — 11 m 215 — 53m 1 4 

- a ^"T6Ll6 + 6 — 4m + m 2 [6 — 4m + m 2 ] 2 30 — 4m + m 2 J m ' 

a 2 — a_ 2 3 T243 , 175 + 563m , Q 89 — 32m 361 — 10m 



48: 



- 361 — 10m -j 4 
2T2+ 30 — 4m-l-m 2 J m ' 



a 32L 8 '6 — 4m-fm 2 [6 — 4m + m 2 ] 2 ' 30 — 4m + m 2 

The evident objection to this form for the coefficients is that it makes the 
several terms very large, and of signs such that they nearly neutralize each 
other, the sum being very much smaller than any of the component terms. 
However it may be possible to remedy this imperfection by admitting three 
terms into the numerators, but, in this way, the problem is indeterminate, 
infinite variety being possible. 

It is remarkable that none of our system of divisors can vanish for any 
real value of m, since the quadratic equations, obtained by equating them to 
zero, have all imaginary roots. In this they differ from the binomial divisors 
met with when the integration is effected in approximations arranged 
according to ascending powers of the disturbing force. 

It is well known that the infinite series, obtained from the development, 
in ascending powers of m, of any fraction whose numerator is an entire 
function of m, and its denominator any integral power of a divisor of the 
previously mentioned series, is convergent, provided that m lies between the 
two square roots of the absolute term of the divisor. Hence any finite 
expression in m, involving these divisors, can be developed in such a series, 
provided that the numerical value of this parameter is less than V 6 . The 
same, however, cannot be asserted when the expression really forms an infinite 

series, as it is in the equation just given for the value of — . Yet, on account 

a„ 

of the simplicity with which these quantities can be expressed in this form, 
a t and a_ x containing each a single term, with an error of the sixth order 
only, this limit is worthy of attention. 

If the parameter m, hitherto employed by the lunar theorists, is taken as 
the quantity in powers of which to expand the value of a* , we shall have 
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m = . And, substituting this value, the principal divisor 6 — 4m + m 

1 — m 

becomes 6 — ■ 16m + 11m 2 . Thus the limits, between which m must be con- 
tained, in order that convergent series may be obtained where this divisor 

intervenes, are ± \*r • When we consider how little, in the case of our 

moon, m exceeds m, it will be plain that the series, in terms of m, are likely 
to be much more convergent than those in terms of m. 



If we inquire what function of m, of the form 

M 



m 



1 -\-aia 



, the quantity 



[6 — 4m + m 2 f 
can be expanded in powers of, with the greatest convergency, it is easily 



found that a should be 



Then putting 



m = 



m 



1 + - r .» 
the divisor 6 — 4m + ni 2 is changed into 

6 + ~«t 2 , 

and there is introduced the additional divisor 1 + — ttt . Here the series will 

o 

be convergent provided ttt is less than 3. It is true the terms involving the 

succeeding divisors 30 — 4m + m 2 , &c, are not benefited by this change of 

parameter, but as they play an inferior role in this matter, I have chosen nt 

as the parameter for the developments of the coefficients a^ in series of 

ascending powers. 

To illustrate this matter, we have, in terms of the parameter ttt, and 

with errors of the sixth order, 



a l "T a — l 



"2 + |tn 



1 + i nt 2 1 + ~ m j 



ttt" 



8j i 



' 5 + T*tt 
1 + 5*' 



27 



1 + 4 



+ 



ttt 



[l + i«t]'' 



ttr 



142 Hill, Researches in the Lunar Theory. 

Expanding these expressions in powers of ttt, we get 

^ = — L W + 2 m ~l> m + 36 m +••••]' 

a i a -i n ttt 2 J_ 5 ttt 3 J- 5 ttt 4 1] ttt 5 _1_ 

Let these series be compared with those which correspond to them in the 
lunar theories of Plana or Delaunay, viz : 

o f 19 a i 131 4 « 383 k i 

m? + - m 3 + -jg w 4 + -^ m 5 + . . . . , 

11 2 i 59 3 , 893 4 , 2855 s , 

The superiority of the former, in convergence and simplicity of numerical 
coefficients, is manifest. 

Much more might be said relative to possible modes of developing the 
coefficients a* in series, but we content ourselves with giving their values 
expanded in powers of m, the series being carried to terms of the ninth order 
inclusive. The denominators of the numerical fractions are written as pro- 
ducts of their prime factors, as, in this form, they can be more readily used, 
the principal labor in performing operations on these series being the reduc- 
tion of the several fractional coefficients, to be added together, to a common 
denominator. 



»i _ 3 2 , 1 ' 8 , 7 4 , 11 5 30749 a 1010521 7 18445871 8 

— — gi m "T" -H- m + Yi~Z m + YiJi m zrmT m -zm^ m — t^r^nr; m 

a 



2 .3 ' 2 2 .3 2 2 12 .3 3 2".3 4 .5 2 10 .3 5 .5 2 



_ 2114557853 9 

2 12 .3«.5 3 ' * ' 

-l-_l? ra 2 -— ~m 3 43 m *_ U m 5 -H m «±» m '4.™ 9 m ! 

2* 1U 3 m 2 2 .3 2 U1 S 3 m W^Ji m + 211.3^.5 m + 2<>.3 3 .5 2 m 



a 



. 13620153029 _o 

4- m 

' 2i 2 .3».5 3 



&2 _ 25 4 , 803 6 , 6109 , _,_ 897599 7 . 237203647 8 _ 44461 407673 , 

a — 2» ' 2'. 3.5 ' 2s.3 2 .5 2 ~r 2 3 .3 3 .5 3 ' 2i«.3 2 .5* '2i5.3*.55.7 

a -2 _ 0rn 4 , 23 5 , 299 6 , 56339 7 , 79400351 8 , 8085846833 9 

% _ 833 6 , 27943 7 , 12275527 8 , 27409853579 9 

9 212.3 "T"2 Il .5.7 "t" 2i°.3 2 .5 2 .7 2 "^ 2 12 .3*.5 3 .7 3 ' ' " 
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— — — — m 6 -I- m 7 -I- m 8 4- — m" 

fl 2«.3 "•" 2'. 3.5 ' 2 8 .3 2 .5 2 .7 ' 2».3*.5 S .7 2 



a_s _ _i_ 

a _ 2«.E 

^4 _3K7 . , 111809667 9 
— 2 16 ~T 2i».3 2 .5.7 2 

a_ 4 _ _28_ 8 , 1576553 9 
<, — 2H.3 ~ t " 2i».3 2 .7 2 ' 



These values being substituted in the equations 

r cos v = Xi • a^ cos 2#r , 
r sin -y == ^ . a^ sin 2ir , 
and the parameter changed to m, we get 

, 25239037 8 7 32931 9 1 ~ 
+ 14929920 m _ 87824800 m * * *J C0S M 

4- ["-- tit 4 4- -- tit 3 4- — tit 6 — 1 8 1L m 7 

~ |_256 m ~ 960 m ~ 28800 m 216000 ,U 

_ 164645363 . _ 11321875589 9 1 . 

552960000 m 19353600000 ul ' ' 'J 

_L T 2 " m e _l_ 80193 m 7 4- 379549 ^v8 j_ 181908179 - "I ft 

+ |_4096 m + 107520 m +100352(r m + 1580544000™ ' ' J C0S 6t 

, T 11347 _ 8 , 2350381 « ~| . ) 

+ Ll966T8 W +"903i680 m ° • • -J COS 8r + . . . j , 

r sin i, — a I ["" tti 2 4- 5 m 3 4- 5 tit 4 — u tit 5 — lou23 m 6 _ 512239 m 7 

r Sm V — a ° | L 8" m + 4" m + 72 m 36 m 331776 W 276480 m 

_ 269023019 8 1 51872119 9 "I . „ 

74649600 ™ 93312000 "* ' * ' ' J 

4. ["-* tit 4 4- — tit 5 4- -— m 6 — i™i tit 7 

~ |_256 ' 480 m ~ 28800 m 432000 '" 

_ 350028 7 8 _ 43885512859 9 "I . . 
H520000 m 58060800000 m ' ' ' ' J 

4- T-^ 9 - tit 6 4- - 4 - 1 m 7 4- ii 19 ! 4 ! tit 8 4- 39 8314169 nt 9 1 sin frr 
' L 12 288 ,H ^ 107520 m ~ r 15052800 ,U "•" 4741632000 m ' ' 'J 

_i_ T 9875 tit 8 4- 32608451 tn 9 1 sin 8-r 4- I 

+ [mm m + 144506880 m j sm or + . . . I . 

Our final differential equations are capable of furnishing only the ratios 
of the coefficients a*, hence we must have recourse to one of the original 
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equations if we wish to determine a as a function of n and (i. By substi- 
tuting the values 

in the differential equation 

\& + 2mD + |-m 2 - * ] « + |m^ = 0, 

L 2 (tts)sJ 2 

we obtain 

Considering only the term of this, for which i = 0, and supposing that the 



&lu 



coefficient of £ in the expansion of ° is denoted by J", we shall have 

{ IIS) ~2 

4j=l + 2m + -|m 2 + -|m 2a ^. 
ao 2 2 a 

For brevity call the right member of this H; then, since 

* = 7-^2 = 4 (1 + m ) 2 . 

(to — n'y w 
we shall have 



„_ [5t] »pfl+-q». ^ 

The value of .ff is readily obtained from the value of -=i given above, and J 

a 

must be found by substituting the values 

in --5— and taking the coefficient of £. We get 

{us)i 

t _ i i f a i + a_ !"| 2 r3 45 t- a t + a _ri 2 15 a t a_! 15 a 2 + a_ 2 ~l 

"*" L a J L4 ^64L a J "^ 8 a 2" a J 

i ^2 i &— 2 I " 85 T <^— 2 1 fi a^a_ j 1 1 e 811 1 a — 1 &1&2 ~r <^— i^— 2 
a L4 a a 2 , J a a 2 . 

■ o *Mjt— 1 1 4 k aia_ 1 1 q a, 2 a — 2 
a a a 

where the terms neglected are, at lowest, of the tenth order with respect to m. 
And, explicitly in terms of this parameter, 

T _ 1 , 21 4 _ 81 6 _ 53 6 _ 2707 - 4201213 8 , 14374939 « 
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By means of which there is obtained 



^=L^JL 1 -t«' i +1- i +S»'-S-'- 


45293 . 
41472 


8761 7 4967441 8 , 14829273 9 
6912 7962624 ' 39813120 


...], 


or, in terms of the parameter m, 




a - = [-S J [ 1 -7» ! + ^*-£> , -a"£<"-- 


79859 , 
124416 m 



4- *®±. m 7 4- i®M96 8 , ii914_o_6_5l_ c 
' 10368 ,U ' 71663616 m "r" 1074954240 ,w 



....]. 



The quantity -^ is usually designated a by the lunar-theorists ; and, 

to make this appear as a factor of the expressions for r cos v and r sin v, it 
would be necessary to multiply all the coefficients by the second factor of the 
preceding expression for a . It seems simpler however to retain a as the 
factor of linear magnitude ; for the astronomers have preferred to derive the 
constant of lunar parallax from direct observation of the moon, or, in other 
words, they have preferred to consider p as a seventh element of the orbit ; 
with this view of the matter, there is no incongruity in making a everywhere 
replace p. 

The expression for a can be obtained in several other ways, which lead 
to more symmetrical formulae, and which also serve for verification of all the 
preceding developments. If, in the preceding equation giving the value of 

-— — in terms of £, we attribute to t the value 0, or, which is equivalent, 

make £ = 1, we shall have u = s = % t . a^ , and, consequently 

= 2 i .[(2;+l + m) 2 + 2m 2 ]a i . 



OS* • a J 

And thus, mindful of the value of % given above, we get 
a = |>-|*r (1 + m) 2 

~ L» 2 J 



X t . \(2i + 1 + m) 2 + 2m 2 ] -« . h t . -/f 



Again the differential equation 



g + ang + ^ = o 



37 
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gives 

-- . y = S t . [(2i + 1 + m) 2 — m 2 ] a, sin (2t + 1) * , 

and, attributing to t the special value — , 

= X i .(- ly (2» + 1) (2» + 1 + m) a, . 



pM-.l)«aJ» 
Whence 



-=te] 



* r 



(1 + m) s 



; { . (_ 1)* (a + 1) (2i + 1 + m) * . fa. (- 1)* * Y 

When j = in the first equation of condition for determining the coeffi- 
cients fy, we get a formula expressing Cin terms of these quantities, viz., 

- %. [(2* + 1 + 2m) 2 + I m 2 ] a 2 + j m 2 2,. a^a.^ , 
or neglecting terms of the eighth and higher orders, 

C= a 2 [~1 + 4m + £ m 2 + (9 + 12m+ 1 m 2 ) -J + (1— 4m + 1 m 2 ) -=± + 9m 2 ±=l\ 

L l l K 2 a 2 , a J 

2 T-i i ( i 9 2 H47 _ 4 1399 * 2047 6 , 3737 7 ~1 

= a 2 [1 + 4m + T m 2 — _- m 4 — -- m° — m + ^i ™ 7 J • 

But the C of Chap. I is obtained by multiplying this C by — v 2 = — -. . 

* _ J f j & j 2 2 (1 + m) 2 

Hence, substituting for a its value, we have 

fc>).[i+».-|«'-»'-s n '-««--Srf-s«'], 

as there stated. 

We propose now to reduce the preceding formulae to numerical results. 
For this purpose we assume 

n = 17325594".06085, 
rt = 1295977".41516, 



which give 



m = — ^— = 0.08084 89338 08312, 

n — n' 

m 2 = 0.00653 65500 97941, 

m 3 = 0.00052 84731 06203, 

m 4 = 0.00004 27264 87183, 

m = 0.08308 81293 65. 
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The numerical value of tit being substituted in the series, we obtain 

a = 0.99909 31419 62 [41* , 

rcosu = ao[l — 0.00718 00394 55 cos 2t 
+ 0.00000 60424 59 cos 4? 
+ 0.00000 00325 76 cos 6* 
+ 0.00000 00001 80 cos 8*] , 

r sin v = a [ 0.01021 14543 96 sin 2* 
+ 0.00000 57148 79 sin 4* 
+ 0.00000 00274 99 ain\r 
+ 0.00000 00001 57 sin 8*]. 

(To be Continued.) 



